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A microscopic derivation of the classical Generalized Constant Coupling (GCC) model for ge- 
ometrically frustrated magnets is presented. Starting from the classical Heisenberg Hamiltonian, 
the partition function for clusters with p = 2, 3, 4, . . . spins in the presence of the inhomogeneous 
symmetry breaking fields (SBF) created by spins outside the unit is calculated. The effective fields 
characterizing the interaction between units naturally arise as averages over the SBF. These effective 
fields are fixed by a self-consistency condition. In the paramagnetic regime, we recover all the results 
previously obtained in a more phenomenological way, which were shown to be in excellent agree- 
ment with Monte Carlo calculations for these lattices. However, this microscopic formulation allows 
us to also study the behavior of the system in the absence of applied magnetic field. It is found that, 
for antiferromagnetic interactions, the equilibrium configuration is a non-collinear configuration in 
which the total magnetization of the unit is zero, and the condition under which such an ordered state 
occurs is also obtained from the calculation. However, frustration inhibits the formation of such an 
state, and the system remains paramagnetic down to K, if only nearest neighbor interactions are 
taken into account, for all the systems considered, in agreement with the now generally accepted 
idea. For completeness, the ferromagnetic case is also studied. Furthermore, as the present method 
is formulated in the real space, it is very easy to study the effect of further perturbations in the Hamil- 
tonian by modifying the form of the SBF. We do this to include the effect of next nearest neighbor 
interactions (NNN) and site dilution by non-magnetic impurities for the pyrochlore lattice. It is found 
that NNN interactions can stabilize a non-collinear ordered state, or ferromagnetic one, depending on 
the relation between NN and NNN interactions, in agreement with mean field calculations, and the 
phase diagram is calculated. However, site dilution is not enough by itself to form such an ordered 
state. 

PACS numbers: 75.10.Hk, 75.25,+z, 75.40.Cx 



I. INTRODUCTION 

The study of the magnetic ordering in geometrically frustrated antiferromagnets (GFAF) is not a 
new problem in Condensed Matter Physics. It was initiated in the fifties by Wannier, Houtappel, and 
AndersonQffl A further milestone was set by VillairE, who first argued that the pyrochlore Heisenberg an- 
tiferromagnet remains paramagnetic down to K. However, there has been a renewed interest during recent 
years on these lattices, mainly because, experimentally, a rich phenomenology has been found at low tem- 
peratures and novel phases have been identified, even though these systems are relatively simple crystalsi@fl. 
In these materials, the elementary unit of the magnetic structure is the triangle, which makes it impossible 
to satisfy all the antiferromagnetic bonds at the same time, with the result of a macroscopically degenerate 
ground state. Among the lattices with the highest degree of frustration are the pyrochlore and kagome ones. 
In the former, the magnetic ions occupy the corners of a 3D arrangement of corner sharing tetrahedra; in the 
latter, the magnetic ions occupy the corners of a 2D arrangement of corner sharing triangles (see Fig. [[]). To 
summarize a very long story, it is found that the susceptibility in materials belonging to this class exhibits 
a high temperature phase in which it follows the Curie-Weiss (CW) law. Below the Curie-Weiss temper- 
ature, geometrical frustration inhibits the formation of a long range ordered (LRO) state, and the system 
remains paramagnetic, even though there are strong correlations between units. This phase is universally 
present in these systems, and it is called the cooperative paramagnetic phase. Finally, at a certain temper- 
ature, Tj, which depends on the particular material and, usually, is well below the CW temperature, there 
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appear non universal phases: some of the systems remain paramagnetic down to the lowest temperature 
reached experimentally!, some of them exhibit non collinear ordered statesBS, or even some of them form 
a spin glass statellllOO, even though the amount of disorder in the structure is very small. It is thus easy to 
understand the amount of theoretical attention these systems have generated0lll001ilJili!30iiil0. 

In the simplest theoretical description of GFAF, the spins on the lattice are regarded as Heisenberg spins 
with only nearest neighbor (NN) interactions. In this picture, it is predicted that the non trivial degeneracy 
of the ground state inhibits the formation of a LRO state, and the system remains paramagnetic down to 
zero temperatureBlll. However, due to the presence of frustration, the NN exchange does not fix an energy 
scale on the problem, and any small perturbation can break the non trivial degeneracy of the ground state 
and lead to some kind of ordered state. Therefore, it is especially important to incorporate these possible 
perturbations in any model that tries to explain the low temperature phases of these systems. Examples of 
perturbations present in real systems are next nearest neighbor (NNN) interactions, small anisotropics, long 
range dipole-dipole interaction, or dilution by non magnetic impurities, to cite some. 




(a) Kagome lattice 



(b) Pyrochlore lattice 
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(c) Linear chain 



Figure 1: The lattices we will consider in this work. A, B, C, . . . denote the sublattices in which the whole lattice is 
subdivided. All these lattices share the following property: if we consider a spin inside a cluster formed by p spins, 
there are p—1 NN spins outside the cluster. Notice also that one spin belonging to one sublattice does not have NN 
belonging to the same sublattice. 

In a series of recent paperscj the present authors have generalized the well known Constant Cou- 
pling method of magnetisml*!!], to deal with geometrically frustrated lattices, in both the classical and 
quantum cases, in the cooperative paramagnetic region. Also, the connection between this technique and the 
so called Phenomenological Renormalization Group (RG) methods for Ising spinsE3Ell0 has been furnished 



in Ref . 33 , where it was shown that the Generalized Constant Coupling (GCC) method can be regarded as 



the zero-th order approximation of one of these RG techniques, for vectorial or scalar spin systems, the so 
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called, Mean Field Renormalization Group method, when correlations outside a cluster are neglected, and 
the anomalous critical exponent of the system is taken to be zero. Therefore, it is not so surprising that 
magnetic properties calculated within the framework of the classical GCC method were found to be in ex- 
cellent agreement with Monte Carlo calculations for those same quantities9EH!, for both the kagome and 
pyrochlore lattices. The main idea behind this technique is very simple, and it is based in the experimental 
and numerical observation that correlations in GFAF are always short ranged§§0E10§l. Therefore, in or- 
der to calculate any thermodynamic quantity, we can start by calculating the partition function of an isolated 
cluster with p spins and, later, add the interactions with the surrounding p — 1 spins outside the cluster (see 
Fig. [j} by means of an effective field which is fixed by a self-consistency condition. This, almost phe- 
nomenological, formulation of the method, allows one to easily compute thermodynamic quantities, mainly 
the susceptibility, in the paramagnetic region, by using a single effective field to characterize all the inter- 
actions with spins outside the cluster. However, it is not very useful to study the critical properties of these 
systems, except for ferromagnetic interactions, which is not as interesting as the antiferromagnetic case, due 
to the absence of geometrical frustration (for Heisenberg spins). The reason is very simple. It is now gen- 
erally accepted that the ordered state of these GFAF, if any, is not a simple antiferromagnetic configuration, 
as in bipartite lattices, but a non-collinear (NC) one, in which the total spin momentum of the unit is iden- 
tically zeroBB Obviously, we cannot characterize such an ordered state by using a single order parameter 
fixed by a single effective field (EF). What is even worse, even we know that we have to introduce different 
order parameters for each of the spins in the cluster, and different (EF) for the neighboring spins outside 
the cluster, it is not easy to intuitively see how these EF enter into the expressions of the order parameters. 
Therefore, we need to go beyond the phenomenological formulation of the GCC method introduced in the 
aforementioned works. This is precisely the intention of this work. We construct the GCC method from 
first principles, by using rigorous identities in terms of averages over finite size clusters. In this way, the EF 
naturally emerge as averages over the microscopic spin variables. This microscopic formulation of the GCC 
method not only allows us to recover all the results obtained Refs. 25 and 27 in the paramagnetic regime 
but, also, gives us the possible ordered configurations of the system, and the conditions for a transition to 
one of these states to occur. Moreover, as the GCC method is constructed in the real space, it is especially 
easy to include further perturbations in the system. Particularly, we study the effect of next nearest neighbor 
(NNN) interactions and site dilution by non-magnetic impurities. 

Apart from the microscopic formulation of the GCC method itself, we present various interesting results 
in this work. We find the possible ordered configurations for systems formed by corner sharing units. These 
systems have the interesting property that the corresponding lattices are formed by units with p = d + 1 
spins (d the dimensionality), and each of the spins in the unit has p — 1 NN outside the unit. Thus, for d = 3, 
we have the pyrochlore lattice; for d = 2, the kagome lattice, and for d = 1, the linear chain. We have 
studied this last case for the sake of completeness, even though this system is not frustrated. It is found that, 
apart from the trivial ferromagnetic configuration for ferromagnetic interactions, the ordered state would 
correspond to a non-collinear configuration given by the condition that the total magnetization of the unit 
is zero (for the especial case of the linear chain this corresponds to standard antiferromagnetic order). For 
NN interactions only, there is no transition to such arj_onlered state at any finite temperature, for none of the 
systems considered, as it is now generally acceptedotS. However, when NNN interactions are taken into 
account, there is a transition to such a NC ordered state for some ranges of the coupling parameters. We 
have only studied the pyrochlore lattice in this case, as we do not think the results of the present method 
for the critical behavior of these systems are reliable in dimension lower than 3 (actually, the results of any 
model based on the idea of finite size clusters), as the long wavelength fluctuations around the NC ordered 
state are not properly taken into account, and they are very likely to be specially important in 2 dimensions 
(kagome lattice). Also, the effect of site dilution is considered in the present work. It is found that, in this 
framework, dilution by itself is not enough to stabilize a NC ordered state. Finally, we also consider both 
dilution and NNN together, and the phase diagrams for different amounts of dilution are calculated. 

The remaining of the paper is organized as follows: In the following Section, we introduce the systems 
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studied in this work, the Hamiltonian of these systems, the notation, and the self-consistency condition that 
fixes the EF. In Section III we present the calculation of the order parameters for finite cluster of size 1 and 
p (p = 4, 3, and 2 for the pyrochlore, kagome, and linear chain, respectively). The connection with the 
original formulation of the GCC method is furnished in Section [V. In Section |V] we study the conditions 
for a transition to a long range ordered state. In Section [v| we include NNN interactions and study how 
this additional interaction modifies the conclusions of Section |v[ In Section VII we develop the formalism 
to deal with site dilution by non-magnetic impurities, and apply it to the present problem. In Section VIII 
we consider simultaneously the effect of both NNN interactions and site dilution. Section IX is devoted to 
the discussion of the results and the conclusions. Finally, we have included two appendixes to clarify some 
of the most involved details of the evaluation of the partition function of the different clusters subject to 
inhomogeneous EF. 



II. PRELIMINARIES 



As stated above, the original formulation of the GCC methodtrn£3 is not well suited to study the critical 
behavior of GFAF Indeed, it is now generally accepted that the ground state of the pyrochlore and kagome 
lattices is fixed by the condition that the total spin momentum of the tetrahedral or triangular unit, respec- 
tively, is zero, though frustration inhibits the formation of such a state, and the system remains paramagnetic 
in the whole temperature range, if only NN interactions are present. Obviously, we cannot characterize such 



an ordered state with a single order parameter, as we did in Refs. [25| and which restricts the applicability 
of the GCC method to the study of the paramagnetic region or, at most, the ferromagnetic case, which is of 
less interest, due to the absence of geometrical frustration. A relatively simple way of circumventing this 
shortcoming of the method has been pointed out in Refs. |4| and 33, and essentially consists on introduc- 
ing different sublattices characterized by different order parameters. The minimum number of sublattices 
we need to introduce is given by the number of spins in the cluster, that is, 3 and 4 for the kagome and 
pyrochlore lattices, respectively (see Fig. [[]). With this subdivision, spins belonging to one sublattice only 
interact with spins belonging to different sublattices. Therefore, the dimensionless Heisenberg Hamiltonian 
of a lattice formed by corner sharing clusters with p spins, and only NN interactions, in the presence of a 
uniform magnetic field can be put in the form 

H = K ^Sia ■ Sjp + H ■ ^Sia, (1) 

where Si a are the classical Heisenberg spins of unit length. The a index labels the sublattice to which the 
considered spin belongs (and takes the values a = A, B, C, . . .), whereas the i index labels the spins belong- 
ing to a given sublattice. K = ^ (K > for ferromagnetic interactions and K < for antiferromagnetic 

ones) and Hq = 4^ , with ho the applied magnetic field. 

The central idea of the GCC method is to calculate the order parameter (for each sublattice) for two 
clusters of different sizes, p and p 1 , replacing the effect of the spins not included in the cluster by fixed 
effective fields, which act as a symmetry breaking field (SBF) and, by using the self-consistency condition 

m p ,(K)=m p (K), (2) 

obtain an equation to determine these effective fields and. thus, the magnetization. The order parameter is 
calculated by making use of the Callen-Suzuki identityEJ E^ 

m ^ = M = \ Tr p en P ) n (3) 

where the partial trace is taken over the set of p variables specified by the finite size cluster Hamiltonian TL P ; 
rh a is the corresponding order parameter, and (. . . } w indicates the usual canonical thermal average over 
the ensemble defined by the complete Hamiltonian Ti. 
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III. EVALUATION OF THE ORDER PARAMETERS IN THE PRESENCE OF THE SBF 

A. 1-spin cluster order parameter 

Let us consider the cluster denned by the spin 1 in sublattice a (see Fig. |2|). The Hamiltonian of this 
cluster is given by 

Hla = Slot \K ^ ^2 % + ^0 = Si a - ii a , (4) 
\ (3^a i J 

where the subindex la makes reference to the fact that the Hamiltonian corresponds to a cluster of 1 spin 
belonging to sublattice a, and £i Q stands for the symmetry breaking field acting on the spin of the 1-spin 
cluster, which belongs to sublattice a. The sum over the i index is performed over the NN of the la spin, 
with z the number of NN in each of the sublattices (z = 2 for both the kagome and pyrochlore lattices). 




(c) Linear chain 

Figure 2: 1-spin clusters for the different sublattices considered in this work. The interactions with spins outside 
the cluster, i.e., those creating the EF, have been represented by dashed bonds. Notice that there are 2 NN on each 
sublattice, different from the one to which the considered spin belongs, outside the cluster. 

By using the Callen-Suzuki identity we can calculate the order parameter in sublattice a for a cluster 
formed by 1 spin 

. / Tr la s la e n ^ \ 
mi " = (Sl " ) = \ Tr la e^ ) H (5) 

As commented above, Tr\ a represents the partial trace with respect to the variables of the spin la, and 
the symbol (. . . } n stands for the usual canonical thermal average taken over the ensemble defined by the 
complete Hamiltonian TL. This equation can be put in the more convenient form 

mia = (V|- ia lnTr la e Hl ^ n . (6) 
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For classical vectorial spins of fixed length, the trace operation is given by 

Z la (K) = Tr la e gla ^ la = [ dQe g ^ a , (7) 



where the integration is performed over the angular degrees of freedom of the spin. This integral is trivially 
evaluated to give 

Z la (K) = 2^ (8) 

and, taking into account the well known relation Vpr = n, where n is a unitary vector in the direction of f, 
we can easily calculate the order parameter for the 1-spin cluster as 

m la = (-L(fia)) n la , (9) 

where L{x) is the Langevin functionH and n\ a is an unitary vector in the direction of the corresponding 
SBF. 

Near a critical point, as we expect the SBF to be very small, we can expand the previous expression and 
keep only linear terms in the SBF 



mia - ^y^- (10) 

By defining the effective fields h' a = J = J (sja}^ ^ 1), we arrive to the final expression for 

the order parameter 

2 v ^ ho 

where we have taken into account that each spin has 2 NN on each sublattice different from the considered 
one. For example, for the pyrochlore lattice with no applied magnetic field we have 

2 - 

rhi a = —(tip + h' 7 + ti s ), (12) 
where a, (3, 7, and 5 can take the values A, B, C, or D. For the kagome lattice 
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where a, (3, and 7 can take the values A, B, or C. For a cluster with 2 spins, we have a linear chain, for 
which 

2 -> 

mi a = — hp, (13) 

where a and (3 can take the values A or B. 

Another interesting point is that if we keep all the terms in the expansion of the order parameter, and 
make use of the decoupling approximation 

{SictSjfi • ■ ■ ^^7)7-^ — {Siajj-i ■ • • (^^7)7-^ j (14) 
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and further assume h' a = rh\ a , we recover the usual self consistent Curie-Weiss equation for the order 
parameter 

m la = L I 2Kj2mi/3 + H I ■ (15) 
\ P¥*» J 

However, this expression would predict a transition to an ordered state (for antiferromagnetic interactions) 
at a finite temperature given by Tjy = — &cw = 2 ^ P ~^ J ^ , which it is well known to be incorrect for 
GFAF. That is one of the reasons why constructing the effective fields by the self-consistency condition (^) 
provides a much better theory for GFAF, as we will see below. 



B. p-spin cluster order parameter 

The calculation of the order parameter for a cluster with p spins in the presence of the SBF created by the 
neighboring spins is a highly non-trivial problem, as we will see in the following. The size of the clusters 
is chosen so that we consider the smallest number of spins that include all the different sublattices. That 
is, p = 3 for a kagome lattice (a triangular cluster), p = 4 for a pyrochlore (a tetrahedral cluster), etc (see 
Fig. ||). The Hamiltonian for such a cluster can be put as 

Hp = K si a • sip + ^s la ■ i pa . (16) 

In this Hamiltonian, the first term represents the interaction between the spins inside the cluster of size p, 
whereas the second term represents the interaction of those same spins with the symmetry breaking fields 

z-1 

f pQ = K^^% + i? , (17) 

P^a i 

created by the p — 1 spins outside the cluster. The sum over i goes over all the NN spins except the ones 
already included in the cluster. 

The order parameter for each of the sublattices included in the cluster, obtained from the p-spin cluster 
can be now put as 



/ TriA,iB,ic...si a e n P 
\ 7>iA,iB,ic...e w P 



(18) 

H 



or, more conveniently 



V? mTr 1AilBiia ..eM . (19) 



In these expressions, TriA,iB,iC... represents the trace over the spin degrees of freedom in thep-spin cluster. 
The reduced trace operation for Heisenberg spins is given by 

Z P {K) = f Y[dn a e n * (20) 
It is impossible to find a closed expression for this partition function of the cluster in the presence of inho- 



mogeneous SBF (the case with no SBF has been studied by Moessner and Berlinsky in Ref. 15). However, 
as we will be mainly interested in the region near the critical point (if any) in the absence of applied field, 
where the SBF are expected to be small, and in the paramagnetic region, where the SBF will be proportional 
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(a) Kagome lattice 



(b) Pyrochlore lattice 



B 



A 



B 



A 

4> 



(c) Linear chain 

Figure 3: p-spin clusters for the lattices considered in this work. Again, we have represented the interactions with the 
spins creating the EF by dashed bonds. Notice that there is 1 spin outside the cluster on each sublattice different from 
the considered spin of the cluster. 



to the applied field and, again, small, it is enough for our purposes to calculate reduced partition function up 
to second order in the SBF. The main steps of the calculation are presented in Appendix |a|, and the result is 
given by 



for p > 1, where 



and 



Z p (K) ~ z° p {K) + 



zl{K) 



6 ^ 



z°JK) 



aqq e2K ' 



o 



zl(K) 



2 -2— 

•q 



smq 



p-2 



cos q sin q 



(21) 



(22) 



(23) 



Expressions for these functions are given in Appendix [B[ Incidentally, the first term in ( |2T| ) corresponds 
to the partition function for non-interacting clusters with no applied field, which were first calculated by 
Moessner and BerlinskyEl. The additional terms in ( |2T| ) allow us to study not only the paramagnetic regime, 
but also the behavior of the system near a critical point, as we will see below. 
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Once we have this reduced partition function we can calculate the order parameter up to first order in 



the SBF by taking into account (19) 



in 



pa 



3 3z°(K) ^ W • 



(24) 



Taking into account the definition of the SBF, and the fact that J (sj Q ) w = J (s)a) w = h' a , after a little 
algebra we arrive at 



in 



pa 



^ | [1 - (P - 1K(*Q] ho + (l-(p- 2)a p (K)) jy p -(p- l)a p (K)h' a J 



(25) 



where we have introduced the function 



a p (K) 



4 (k) ■ 



(26) 



IV. CALCULATION OF THE SUSCEPTIBILITY: CONNECTION WITH THE ORIGINAL 
FORMULATION OF THE GCC METHOD 

In the following, we will find more useful to express all the thermodynamic quantities in terms of the 
dimensionless temperature 



1 T 

\K\~U~\ 



In terms of this parameter, we have 



and 



in 



1 



pa 



3\J\T 



1 - (p - l)A p (f )J h + (l - (p - 2)A p (T)j Y,h'p-(j>- l)A p {T)h' a ) , (29) 



(27) 



(28) 



where 

A P (T) = a p (f) (30) 

for ferromagnetic interactions and 

A p {f) = a p {-f) (31) 
for antiferromagnetic ones. The self consistency condition (see eq. (^)) 

rh la = m pa (32) 
leads to the following system of linear equations 

(p - 1) A p (T) h' a +[l + (p- 2) A p (T )] Y, % = ~(P - 1) Ko, (33) 
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that has the solution 



h' a = h' 



1 + (p - l)A p (T) 



ho. 



Substituting back in (g|) we find 



rn 



pa 



1 l-(p-l)A p (T) % 



3\J\T l + (p-l)A p (T) 



that is, the susceptibility is given by 



gcc = 1 _ l-(p-l)A p (T) 

P 3|J|f l + (p-l)A p (T) 
By direct substitution, it is simple to verify the relation 

e p (f) = -(p-l)A p (T), 

where e p (T) is defined as@ 

2T 2 d 



P dT 



lnz°(T)-l 



so we can finally put 



ncc = 1 l + e P C0 



3\J\T l-e p (T) 



in complete agreement with the result in Ref. I 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



V. CALCULATION OF THE CRITICAL POINTS 



In the absence of applied field, the system of equations ([33|) takes the simpler form 

(p - 1) A p (f ) K' a +[l + (p- 2)A p (f)] £ h' p = 0, 



and it has two different non-trivial solutions. The first one occurs if 

A P (f) = 1 

or, equivalently, 

e p (f) = -(p-l), 
and corresponds to an ordered state in which 

rh a = 0. 



(40) 



(41) 



(42) 



(43) 



We will call this kind of state a non-collinear (NC) ordered state. Actually, this condition includes, as 
particular cases, states in which spins in different sublattices are antiferromagnetically aligned with each 
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one. The second solution of (40) occurs if 



other, that is, collinear states. However, the most general state described by this condition is a non-collinear 

(44) 



A P (K) 



1 



p — 1 



or, equivalently, 

e p (f) = 1, 

and corresponds to an ordered state in which 

fh a = rhp = fh^ = . . . 

that is, ferromagnetic order. 



Vex,/?, 7, . 



(45) 



(46) 




Figure 4: Graphical representation of the conditions for (a) Non-collinear order and (b) Ferromagnetic order. 

The previous condition for the existence at any finite temperature of a ferromagnetic ordered state is 
equivalent to studying if the function G P (T) = e p {T) — 1 changes its sign for any finite value of T. This 
function has been depicted in Fig. [4(b) for the pyrochlore lattice, kagome lattice, and linear chain. We can 



see that the present method correctly predicts the existence of a transition for the pyrochlore lattice for a 
value finite value of T {K c = 0.71), and no finite T transition for the linear chain, but incorrectly predicts 
a transition for the kagome lattice (K c = 1.33). This is due to the fact that the long wavelength fluctuations 
around the ferromagnetic state are not taken into account in the present method, and these fluctuations 
destroy the ordered state in 2 dimensions. On the other hand, in the linear chain case, even short range 
correlations, properly described in the model, are enough to destroy any form of long range order. 

The condition for the existence of a non-collinear ordered state is equivalent to study if the function 
F p (T) = e p (T) + (p — 1) changes its sign for any finite value of T. This function has been depicted in 
Fig. |4(a)| . We can see that the model correctly predicts that there is no transition to an ordered state for any 
of the systems considered. 
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For comparison, we have also applied the present GCC method to the Heisenberg model with NN in- 
teractions on a simple cubic (s.c.) lattice. In the framework of the CW theory, it is predicted that both the 
s.c. and the pyrochlore lattices should have the same transition temperature to a LRO state, as they have the 
same number of NN. For the s.c, it is enough, for our purposes, to consider 2 different sublattices. Each 
spin in one sublattice has 6 NN in the other sublattice. Once we take into account this difference with re- 
spect to the pyrochlore lattice (which has p—1 NN, as mentioned above), we obtain that, for ferromagnetic 
NN interactions, there is a transition to a ferromagnetic state for K c = 0.615, to be compared with more 
refined estimates of this quantity from high-precision MC calculations!*], K c = 0.693. For antiferromag- 
netic interactions, a transition to an antiferromagnetic (AF) state occurs at K c = —0.615. Therefore, we 
can see that the present method is able to distinguish between the different topologies of the lattices, even 
when they have the same number of NN. Furthermore, the estimation of the critical point for the s.c. is very 
good (up to -10%), even though we have considered the smallest possible cluster that can describe the AF 
ordered state of the s.c. lattice. 

As mentioned in the Introduction, one consequence of the presence of geometrical frustration in the 
system is that the NN exchange does not define an energy scale in the problem. Therefore, any perturbation 
is a strong perturbation, and can break the non-trivial degeneracy of the ground state, selecting one ordered 
state. In the following sections, we will study if NNN interactions and site dilution by non-magnetic im- 
purities, which are always present in real systems, really break this non-trivial degeneracy, that is, if these 
effects are able to select an ordered state at finite temperature. Furthermore, we will focus the discussion on 
the pyrochlore lattice, for the reasons explained in the Introduction. 



VI. EFFECT OF NNN INTERACTIONS 



In order to include the effect of NNN interactions, we only need to redefine the SBF in the following 



way 



1 



iftl 



3|Ji|T 



^o + 2(l + z' 

P^a 



(47) 



P^a i Pf^a i 

where z' stands for the number of NNN on each sublattice (z' = 2 for the kagome lattice and z' = 4 for 

h, 
■li- 



the, pyrochlore lattice) and A = with J\ the NN exchange coupling and J 2 the NNN exchange coupling. 



For the cluster with p spins 



z-l 



1 



i/3) 



P^a i 



P^a i 



2,\Jl\T 



h + (l + z'X)J2hp 

P^a 



and 



P^a 



(p - l)h + (1 + z' X) 



(p-l)h' a + (p-2)J2h'p 

P+a 



3|Ji|T 

Therefore, the self-consistency condition that fixes the effective fields is given by 

(1 + z' X)(p - l)A p (T)h' a + [l + (1 + z' X)(p - 2)A p (f)] Y,h' p = -ip- l)A p (f)h , 



P^a 



and its solution is given by 



hi 



I + (I + z> X)(p - l)A p (T) 



hn. 



(48) 



(49) 



(50) 



(51) 
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Substituting back into the expression of the 1-spin cluster order parameter we find the expression for the 
susceptibility 



Y gCC 
A,p 



1 



l-(p-l)A p (T) 



3 1 Ji| n + (i + A)(p - i)^4p(r) 



(52) 



or, in terms of the e„ function defined above 



X 



gcc 



1 + e p (T) 



3|Ji|Tl-(l + z' \)e p (T) 



(53) 



This susceptibility has been depicted in Fig. ||] for different values of A and antiferromagnetic NN interac- 
tions. 
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Figure 5: Effect of NNN interactions on the susceptibility of the pyrochlore lattice. 



In zero applied field, the system ( |50P has two non-trivial solutions, similar to the ones found for NN 
interactions only. The ferromagnetic solution occurs if 



and the non-collinear solution occurs if 



>( T ) 



(i + z'xy 



(p-i) 



l + z'X 



(54) 



(55) 



These conditions have been studied and the main results are summarized in Fig. ||. As we can see from 
those figures, NNN interactions can stabilize different ordered states depending on the value of A and the 
sign of the NN interaction. For antiferromagnetic NN interactions and \J^\ < 0.5 | J\\, the system remains 
paramagnetic down to K. For antiferromagnetic NN and NNN interactions, a NC ordered estate is selected 
at finite temperature for A > \. However, for antiferromagnetic NN interactions, and ferromagnetic NNN 
interactions, a ferromagnetic ordered state is selected for A < —\. On the other hand, the case of ferro- 
magnetic NN interactions is slightly more complex. For ferromagnetic NNN interactions, a ferromagnetic 
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(a) Ji < (b) Ji > 



Figure 6: Phase diagram for the pyrochlore lattice with NNN interactions. The curves are obtained for a fixed value 
ofJi. 



ordered state is selected for — ^ < A. In the range — | < A < — g, the system remains paramagnetic down 
to K and, in the range A < — i, a NC ordered state is selected. 

In any case, from the analysis of this section, it seems very unlikely that the noncollinear ordered states 
experimentally found in some pyrochlore systems are due to NNN interactionsQO, as in real pyrochlore 
systems | J2I ~ 0.1 1 J\\, a value well inside the region of the phase diagram in which the system is found to 
remain paramagnetic down to K. 

It is also important to stress that the present treatment of NNN interactions is only approximate. In fact, 
in order to treat this problem consistently, we should include new sublattices such as spins in one sublattice 
do notinteract with spins on the same sublattice. In doing so, it can happen that new kinds of ordered states 
arise&til However, our main intention here is to check if the NC state given by (43) is selected by NNN 
interactions. The general problem with more sublattices will be studied elsewhere. 



VII. EFFECTS OF SITE DILUTION BY NON-MAGNETIC IMPURITIES 



In order to study the effect of site dilution on the susceptibility and on the conditions for the appearance 
of an ordered state, we need to average both the SBF and the order parameters over the distribution of 
non-magnetic impurities. For simplicity, we will assume that the distribution of non-magnetic impurities is 
purely random. Under this condition, the number of units with q spins, for a lattice which, in the absence of 
dilution, is formed by units with p spins, for a concentration of non-magnetic impurities x, is given byi 



pp( x ) = Q (1 - x ) q x p - q . 



(56) 
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However, this distribution is normalized with respect to units, whereas we need a distribution function 
normalized with respect to a single spin. The adequate distribution function, for the SBF, is then 

K(x) = (57) 

and verifies the important relations we will use below 

= 1, (58) 

9=1 



(59) 



9=1 



The corresponding distribution function for averaging over the order parameter, however, has a different 
normalization. Indeed, if we average the SBF for the 1— spin cluster with ferromagnetic interactions with 



respect to (57) we have 



9=1 



hp 
\J\f 



+ 



2(P-1)(1 
I JIT 



-h' 



(60) 



in agreement with the intuitive idea that, upon dilution, the applied field is unaffected, and the EF changes 
according with the reduction of the average coordination number by a factor (1 — x), indicating that the 
distribution function for the SBF is properly normalized. If we now average the order parameter with respect 
to the same distribution function, we find 



x m 



o 



qa 



3\J\T 



+ 



2(p-l)(l^- 
3\J\T 



-h>, 



(61) 



and this result is incorrect, as it does not take into account the reduction in the total number of magnetic 
ions in the lattice upon dilution. Therefore, the correctly normalized distribution function we have to use in 
order to average over the order parameter is given by 



Qg(s) = (1 - xyp»(x). 
In the following, we will use the following notation to denote averages with respect to dilution 



(62) 



(63) 



9=1 



where f q is any quantity, and it will be understood that we are using the distribution ( |57[ ) when averaging 
over SBF and the distribution ( |62| ) when averaging over the order parameter. 

Furthermore, for the sake of clarity, we will consider the antiferromagnetic case with no applied field 
separately from the paramagnetic and ferromagnetic with no applied field cases. In the paramagnetic case 
and ferromagnetic case in the absence of an applied field, we have seen that the system can be described by 
a single effective field, h'. The average with respect to dilution of the SBF for the 1-spin cluster, as we have 
seen above, is given by 



2(p-l)(l 
UIT 



'-h' + 



hp 



(64) 
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where we have only considered the contribution of NN to the SBF. Now, the averaged 1-spin order parameter 
is given by 



1 — x 



Fit 



v 3\J\T 
Analogously, for the p-spin cluster we have 



h + 2(p - - x) h! 



qa 



(p-1) (!-*)_£, , h 



\J\T \J\T 



(65) 



(66) 



^ g-1 
p |J|T 



(p-l)(l-x)K' + K ) , 



and the corresponding averaged order parameter is given by 

r ^ i 1 — X 



3\J\T 



l-Ap(T) /lo + Cp-l)^-^ 7 



where 



I p (T) = ^^(x) (g-l)^(T). 
? =1 



Therefore the averaged self-consistency condition 



leads us to the following equation for the effective field h! 



h + 2(p - 1)(1 - s)# = (1 - A p ) (h + {p- 1)(1 - 



that has the solution 



A, 



(p-l)(l-x)(l + A p ) 
and the corresponding susceptibility is found to be 

1 - x 1 - A, 



ho, 



[xi 



3\J\Tl + A p : 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



which is of the same form as (36), except for that (p — 1)A P is replaced by the corresponding average over 
dilution. Interestingly, it is easily shown that this average can be expressed in terms of an average over the 
e p function introduced before. Indeed, if we define 



2T 2 d ~ 
:,, — — r ^^ln.Zp(T, x) — 1, 



p(l - x) dT 

where z® is the average of z® with respect to dilution, given by 



(74) 



(75) 



9=1 



17 



we have the identity 

and we can put the averaged susceptibility in the form 

1 — x 1 + e Tj 



[xi 



3|J|Tl-e F 



Incidentally, this is precisely the susceptibility calculated by the present authors in Ref. [27 
all the results presented there for the paramagnetic case. 



(76) 



(77) 



so we recover 



i r 



i r 



1.5 - 



Paramagnetic 




i l i 
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(a) Ferromagnetic phase diagram 



(b) Condition for non-collinear order 



Figure 7: (a) Phase diagram for the ferromagnetic pyrochlore with site dilution, (b) Study of the condition for the 
existence of non-collinear order. F4 stands for the function F p = e p — (p — 1)(1 — x) forp = 4 (see text). 



In the absence of applied field, it can be easily shown from eq. (71) that the condition for the formation 
of a ferromagnetic ordered state is given by 

e p = 1. (78) 

To check the consistency of the results, we have studied this condition for different amounts of dilution for 
the ferromagnetic pyrochlore, and the dependence of the transition temperature with x is depicted in Fig. ^. 
We can see that there is a critical value x c above which there is no transition. This value is in reasonable 
agreement with the value of the percolation threshold for the pyrochlore lattice, x c ~ 0.6, whereas our 
model predicts a value x c ~ 0.66. This result is reasonable, as the transition temperatures predicted by our 
model are slightly larger than the real ones, due to neglecting long wavelength fluctuations. 

Let us now consider the antiferromagnetic case with no applied field, and let us calculate the condition 
for the existence of non-collinear order upon dilution. In this case, we know from the analysis in the 
previous sections that the effective fields h' a verify the relation 



-h' 



(79) 
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at the transition point, if any. Therefore, we have the following expressions for the dilution averaged SBF 
in this case 



qa 



2K 
I JIT' 



and 



p+a 



I jit' 



The order parameter for a unit with q < p spins is then given by 

1 



qa 



3\J\T 



+ h' a . 



(80) 



(81) 



(82) 



If we now naively average the self consistency condition as above, we would not be able to obtain the 
condition for the transition to a non-collinear state in terms of the A p or e p functions introduced before, as 
a (q — 1) prefactor is lacking in the expression for m ga . Thus, we construct the following averages instead 



[{q - 1) m la ] p = — — ~ fc a 



and 



[((? - 1 )m 



3|J|T 



((p-l)(l-x)+^ p ) h' a . 



(83) 



(84) 



Equating both quantities, we easily obtain the condition for the existence of non-collinear order upon site 
dilution 



-(p-l)(l 



(85) 



where we have made use of relation (76). 

As can be seen from Fig. 7(b) , the function F p = e p — (p — 1)(1 — x) does not change its sign for any 
value of x, indicating that a non-collinear ordered state is not selected by site dilution, in contrast with the 
qualitative picture suggested in other works. 



VIII. COMBINED EFFECT OF NNN INTERACTIONS AND SITE DILUTION 

Of course, in real systems, the effects analyzed above are always present simultaneously. Therefore, it 
is important to study how the conclusions reached in the previous sections are modified when we simul- 
taneously consider both NNN interactions and site dilution. We will only present the main results of the 
analysis, as the derivation of the final expressions follows very closely the previous derivations. 

The averaged susceptibility is obtained to be 

The condition for the existence of a ferromagnetic ordered state is given by e p = 1 A , whereas the 

corresponding one for a non-collinear ordered state is given by e p = — ■ The phase diagrams, 

for different concentrations of non-magnetic impurities are depicted in Fig. ||. Also, in Fig. ^ we have 
depicted the values of X c (x), that is, the value of ^ above which a non-collinear ordered state is formed at 
a finite temperature, for a given value of x. We have only depicted this quantity for antiferromagnetic NN 
interactions, as the corresponding A c for ferromagnetic NN interactions does not change upon dilution. 
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(a) Ji < (b) Ji > 



Figure 8: Phase diagrams with both NNN interactions and site dilution. 
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Figure 9: Minimum value of A for which a NC ordered state is formed for different concentrations of non-magnetic 
impurities. 

IX. SUMMARY AND CONCLUSIONS 

In this work we have presented a microscopic formulation of the Generalized Constant Coupling method 
for geometrically frustrated magnets. By introducing different Bravais sublattices characterized by different 
order parameters we are able to describe more general kinds of long range ordered states, including non- 
collinear ones. These order parameters are evaluated by making use of the Suzuki-Callen identity for finite 
size clusters, in terms of symmetry breaking fields arising from interactions with units surrounding the 
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considered one. The symmetry breaking fields are easily related with the phenomenological effective fields 
introduced in the original formulation of the generalized constant coupling method, in terms of averages 
over the neighboring spin variables outside the unit, and they are fixed by a self-consistency condition. 
Though we have not stressed this point, this self-consistency condition is, in fact, a particular case of the 
more general scaling hypothesis between the order parameters for finite clusters of different sizesti 

When a magnetic field is applied to the system, this method allows us to easily compute the suscepti- 
bility, and we recover previous results obtained in the phenomenological formulation of the GCC method, 
which were shown to be in excellent agreement with MC results for both the kagome and pyrochlore lat- 
tice. However, the present formulation, allows us to also study the critical behavior of the system when 
no applied magnetic field is present. The possible ordered configurations for these systems naturally arise 
from the equations that fix the order parameters. In completely agreement with the accepted idea, for fer- 
romagnetic interactions a ferromagnetic configuration in which all the spins point along the same direction 
is found. For antiferromagnetic interactions, a non-collinear state in which the total magnetization of each 
cluster is identically zero is found. Also, the conditions for the formation of such ordered states arise from 
the calculation. It is found that, for nearest neighbor antiferromagnetic interactions, frustration inhibits the 
formation of a long range ordered state, and the system remains paramagnetic down to K. We have briefly 
compared the results for these frustrated geometries with the results for a simple cubic lattice obtained with 
this same method, and we found that the generalized constant coupling method correctly distinguishes be- 
tween different topologies of the lattice, even in case where the number of nearest neighbors is the same, 
as it happens with the pyrochlore and cubic lattices. Another advantage of the present formulation is that it 
is constructed in the real space, which makes it very easy to include additional perturbations to the system, 
by modifying the form of the symmetry breaking fields. Particularly, we have studied the effect of next 
nearest neighbor interactions and site dilution by non-magnetic impurities for the pyrochlore lattice. In the 
former case, the phase diagram has been calculated. It is found that, for certain values of the ratio A = y-, 
a non-collinear or ferromagnetic ordered state can be stabilized. However, the minimum value of A above 
which a non-collinear configuration is estable seems to be too large to explain the formation of long range 
order found in some systems, notably in Gd2Ti2O70. Moreover, it seems clear from neutron scattering that 
the ordered configurationin this case is not compatible with the condition that the total magnetization of the 
tetrahedral cluster is zeroll^l. However, the main perturbation in this system seems to come from long range 
dipolar interactions, which are highly anysotropic. Further theoretical work in this direction, in the present 
framework, is still necessary. 

Regarding the effects of site dilution, it has been argued that this kind of effect, always present in real 
materials, could break the non-trivial degeneracy of the ground state and cause a transition to some kind 
of long range ordered state at finite temperature. However, in our calculation, we find that this is not case: 
Site dilution by itself does not induce a transition, for any amount of dilution. It is necessary to include 
simultaneously both dilution and next nearest neighbor interactions. 

A word of caution must be said here, however. In the last part of this work, we have focused our attention 
on studying if a non-collinear ordered state given by the rule that the total magnetization of the tetrahedral 
unit is zero is estabilized by any of these perturbations. However, these are not the xirily possible ordered 
configurations when we include these effects. For example, Reimers and coworkers^ showed that, in the 
framework of mean field theory, next nearest neighbor interactions can induce both non-collinear ordered 
states as the ones studied in this work, or ordered states characterized by an inconmensurate wave vector. 
It is not easy to describe these inconmesurate states with a real space method, and it could be the case that 
site dilution would lead to such an inconmesurate state. Therefore, more theoretical work in that direction 
should be necesary before ruling out the possibility of a long range ordered state induced by small amounts 
of site dilution or, what wouldbe even more interesting, the possibility of a transition to some kind of spin 
glass state induced by dilutionc . 

In conclusion, the formulation of the generalized constant coupling method presented in this work pro- 
vides us a conceptual framework in which both thermodynamic quantities and the critical behavior of geo- 
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metrically frustrated magnets can be properly described. The method can be easily generalized to deal with 
further perturbations always present in real systems. We have presented the classical limit of the method. 
However, in order to compare the calculated quantities with experimental data for these systems, the cor- 
responding quantum generalized constant coupling method should be used. This issue, the occurrence of 
inconmensurate states due to dilution, and the effect of the inclusion of long range dipolar interactions in 
the Hamiltonian are open issues which deserve further theoretical work. 

Appendix A: EVALUATION OF THE P-SPIN CLUSTER PARTITION FUNCTION. 



In this appendix we present the main step for the evaluation of the integral (g0() 

z v(K) = f Hdn a e H ^ (Al) 



which can be put in the more convenient form 



Z P (K)= d 3 S Hd 3 s la 5[S-^2 

si a exp \k y~] s la ■ Sip + 

^ 1 Slot ' Cpa J 

ot \ a J \ a^f3 a J 



(A2) 



6(g) = — 1— / d 3 qe^, (A3) 



By using the Fourier transform representation of the Dirac-J 

1 

W) 

we can further put 

Z P (K) = ^ J d 3 qj d 3 Se% s2 J d 3 s a e s ^, (A4) 

where we have introduced the vector a = £ pa — iq, and made use of the fact 

1 S 2 

with S the total spin of the unit, up to an unimportant global constant in the Hamiltonian. The integral over 
s a can be done as follows. Taking into account 

a ■ s = a x sin 6 cos ip + a y sin sin ip + a z cos 6, (A6) 

where we have assumed that the length of s a is 1 (this can be always done by redefining K — > Ks^ with 
so the length of the spin), the integral over s a can be put as 



J* d6 dp sin 9 e a * cos e+a " sin 9 sin v+a * sin 8 cos * = 2tt d6 sin 9 e a * cos 6 I„ ( + a\ sin I 



/ 1 or I i o\ sinh Va^TV sinh J a 2 x + a 2 + a 2 z { 

2tt / dxe~ ax I [byl - x 2 ) = 4vr / = 4vr / — = 4vr -, (A7) 



V« 2 + b 2 ' . / a 2 + a 2 + a 2 



where a a = q + i £, pa (we have factored out a (-1) in the definition of a a ). 
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Next, we can do the integral over S 



[ SSe^ e* § = ^-^ZLe-^-iO (A8) 
J {-Kf/ 2 



which only converges for K < 0, and we are left with an integral over q 

sin \J q 2 + 2iq ■ t[ pa - £j 



M K )~ / d'qezK [I j = , 

a y q 2 + 2iq '■ £ pa - Q a 



(A9) 



which cannot be calculated in a closed form. An interesting point here is that, if we take £ pa = 0, we 
recover the partition function calculated in the original GCC model. The only thing we can do at this point 
is to expand the partition function in powers of the SBF, up to second order, as the rest of terms will be 
assumed to be small near the critical point once we take the thermal average. With a little algebra it is easy 
to show that 



sin J q 2 + 2iq '• £ pa - & 



" a ~ fo(q) + * fl(q) Z P a cos 6 a + f 2 (q) £ 2 + f 3 (q) f va cos 2 $ a , (A10) 



1 2 + 2iq-£ pa -e 



,pa wo u a -r J'2\HJ ^pa t J6\HJ Spa 



■pa 



where 9 a is the angle defined by £ pa and q, and 



Mq) = = (All) 



cosg sing 

Mq) = o~ ( A12) 

q q 2 



sing _ cos, 
2g 3 2q 2 



, 3 cosg 3 sing sing 

hiq) = ^r-^r + ^q- (A14) 



By making use of this result, we can put 



sin^g2 + 2zg--^-^ a 
1 1 — / = - fo +fo h t pa + fo h 2^ COS 



• \/q 2 + 2iq-^ pa -Q a 



pa 



f P<2 f 2 1 
~ ° 2 1 ^2 COS COS e ^pa£pf3 + i /o _1 /l X] COS ^P"" ( A15 ) 

a^/3 a 

The expressions of the form £ pa cos 9 a are the projection of £ pa over g, and can be put in the more convenient 
form 

Spa COS 6 a = i pa sin 6 cos ip + ^ Q sin sin <£> + £ pa cos 0. (A 16) 
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If we consider the imaginary part of the previous series, and integrate over the angular part of q in (AS), we 
have integrals of the form 



/•27T 

dO / dip sin 9 £ pa cos 6 a , 
(I Jo 



which are easily shown to be zero. 
We also have integrals of the form 



f 2n 

dO I dtp sin 8 (t; pa cos 6 a)' 
(I Jo 



and 



/•2tt 

dO j dip sin £ pa cos 9 a £ P( g cos dp 
o Jo 



4tt 
T L 



(e pa ) 2 + (a y pa ) 2 + a 



pa j 



3 



4tt 



tx tx , ty tV i ez « 
>,paSp(3 ~ ^pa^p/3 ' ^>pa^pf3 



47T 



(A17) 



(A18) 



£ pa • (A19) 



Substituting these results back in ( A15) we are left with the integral over q 



Z p (K) 



a#3 



which, taking into account 



h + 



h fo 



(A21) 



can be simplified to the form 



rp-2 f 2 
JO Jl 



(All) 



Each of the integrals can be put as a closed expressions, and are given in the next section. For the moment, 
we will use the notation 



poo 2 

z° p (K)= / dqq^ehtfV 
Jo 



and 



Finally, up to second order in the SBF, we can put 



d qq 2 e^r - 2 fi 



Z p {K) ~ z»(K) + 



6 



(A23) 



(A24) 



(A25) 
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Appendix B: EXPRESSIONS FOR Zj,(K) AND Z° P {K) 



The Zp(K) function for different values of p and antiferromagnetic interactions (K < 0) is given as 
follows 



z° 1 (K) = e-, 



(Bl) 



K 



(B2) 



3 erf 



4(K) 



+ erf ^3 



(-K) 5 / 2 



(B3) 



2 erf ( V^2K) - erf (V=8K) 3 + e 



SK 



4e 



-2-KT 



(-K) 3 / 2 



8vr K 2 



(B4) 



In these expressions, erf represents the error function. It is important to notice that, if we make the re- 
placement K — > — we immediately recover the expressions of Ref. 15. For ferromagnetic interactions 
(K > 0), the expressions are completely similar to the ones above, with the substitutions \J — K —* yK 



and erf — > erfi, with erfi(x) 



erf (i x) 



the imaginary error function. To quote an example 



2 erfi ( V2K ) - erfi ( V8K 



3 + e 



8K 



4e 



-2K 



K 3 / 2 V8^K 2 
The Zp (K) functions are obtained to be, for antiferromagnetic interactions 



z \ (K) = (sinh K-K cosh K) , 



(B5) 



(B6) 



e 9K f e -AK 



1) 



K 2 



7T l + K 

2 f-iT) 3 / 2 



3 erf 



-if 



erf 3 



-if 



(B7) 



7T 4if + 3 
2 (-X)5/2 



2 erf ( V-2K) - erf ( y/-8K 



ri +K)(3 + e 8K - 4e 2K ) 



(B8) 



The corresponding expressions for the ferromagnetic case can be obtained by using the prescription indi- 
cated above. It is important to notice that Zp {K) is not defined for p = 1 but it easy to see that this function 
does not appear in any calculation, as it corresponds to non-interacting spins. 
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